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Abstract. 

The radial structure of the continuous spectrum of shear Alfven modes is calculated 
in the presence of a magnetic island in tokamak plasmas. Modes with the same helicity 
of the magnetic island are considered in a slab model approximation. In this framework, 
with an appropriate rotation of the coordinates the problem reduces to 2 dimensions. 
Geometrical effects due to the shape of the flux surface's cross section are retained to 
all orders. On the other hand, we keep only curvature effects responsible of the beta- 
induced gap in the low-frequency part of the continuous spectrum. New continuum 
accumulation points are found at the O-point of the magnetic island. The beta-induced 
Alfven Eigenmodes (BAE) continuum accumulation point is found to be positioned at 
the separatrix flux surface. The most remarkable result is the nonlinear modification 
of the BAE continuum accumulation point frequency. 
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1. Introduction 

Plasma stability is one of the crucial issues for fusion devices. Shear Alfven instabilities 
can resonate with energetic particles and are therefore particularly dangerous for a 
burning plasma [H El |3j HJ [5l [6]. Shear Alfven waves (SAW) are electromagnetic 
plasma waves propagating with the characteristic Alfven velocity va = B/y/Ang (B 
is the magnetic field and g the mass density of the plasma). One of the main 
damping mechanisms of shear-Alfven modes in nonuniform plasmas is continuum 
damping [7J El [9], due to singular structures that are formed at the SAW resonant 
surfaces. Due to nonuniformities along the field lines in toroidal geometry, gaps appear 
in the SAW continuous spectrum [10J. The mechanism is similar to that which creates 
forbidden energy bands for an electron traveling in a periodic lattice [HI [121 H3] • Two 
types of collective shear Alfven instabilities exist in tokamak plasmas [14J: Energetic 
Particle continuum Modes (EPM) [15], with frequency determined by fast particle 
characteristic motions, and discrete Alfven Eigenmodes (AE), with frequency inside 
SAW continuum gaps Discrete AE are practically unaffected by continuum 

damping [5l [6]. For this reason, the importance of understanding the continuous 
spectrum structure is clear, if one faces the tokamak stability problem and its potential 
impact on reaching the ignition condition. 

SAW nonlinear dynamics is a research topic with many open issues. In the special 
case of uniform plasmas a peculiar state exists, called the Alfvenic state, which arises 
whenever we can assume ideal MHD, plasma incompressibility and the validity of 
the Walen relation (5v/va = ±5B/B, where 5v and 5B are the perturbed velocity 
and field) [HI [T71 [18]. In such a case, the nonlinear effects due to Maxwell and 
Reynolds stresses cancel and give a self-consistent nonlinear state. On the other hand, in 
nonuniform tokamak plasmas, with non-ideal effects such as resistivity or finite plasma 
compressibility, the Alfvenic state conditions are violated and SAW are characterized 
by a rich nonlinear dynamics. The study of the nonlinear modification of the SAW 
continuous spectrum due to MHD fluctuations can be of practical interest in realistic 
tokamak scenarios. In fact, in a tokamak plasma, the SAW continuous spectrum is 
modified by the interaction with low-frequency MHD fluctuations, such as magnetic 
islands, which are formed when the original sheared equilibrium magnetic field lines 
break due to non ideal effects (in particular finite resistivity) and reconnect with different 
magnetic topology JT9] . Since the typical island frequency and growth rate are much 
lower than the SAW oscillation frequency, we can model the equilibrium magnetic field 
as the sum of a tokamak axisymmetric part plus a quasi-static helical distortion due to 
the magnetic island. 

Here, we derive the fluid theoretical description of the SAW continuum structure 
in the presence of a finite size magnetic island [201 [21] in finite-/? plasmas, keeping 
into account only toroidal effects due to geodesic curvature, which are responsible 
for the beta-induced Alfven Eigenmodes (BAE) gap in the low frequency part of the 
SAW continuous spectrum [221 [23], [21]. We adopt a linear ideal MHD model, and we 
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consider only shear Alfven modes with the same helicity of the magnetic island. In this 
framework, with an appropriate rotation of the coordinates the problem reduces to 2 
dimensions. Since we are interested in the structures of the SAW continuous spectrum, 
characterized by local "radial" singular behaviors, the problem can be further reduced to 
one dimension. The local differential equation describing the nonlinear SAW continuum 
structure is solved numerically with a shooting method code in the whole spatial range of 
interest (both outside and inside the island), and analytically near the O-point and at the 
island separatrix, where the equilibrium quantities can be approximated in simple form. 
Note that what we are calculating is essentially the distortion of the SAW continuous 
spectrum in a one-dimensional system (ID-slab or cylinder) due to the presence of a 
magnetic island. For this reason, for given island helicity, we refer to the structures 
derived here as "2D continuous spectrum of SAW" in the presence of a magnetic island. 
Note also that the problem investigated in this work is somewhat different with respect 
to that studied in [251 . where modes with different helicities from that of the magnetic 
island are considerectjj, yielding 3D structures of the SAW continuous spectrum. 

As suggested by the magnetic field line helicity behavior [261 12Z] , in an equilibrium 
with a magnetic island the separatrix flux surface plays an important role, hosting 
the BAE Continuum Accumulation Point (BAE-CAP), which - without island - was 
positioned at the rational surface. The degeneracy in frequency of even and odd modes 
is removed by inhomogeneities along the field lines, causing a splitting between the 
continuous spectrum branches with different parities. As a consequence, the BAE- 
CAP is also split in different frequencies, depending on the mode number and parity 
of the eigenfunction. Several branches of the nonlinear SAW continuous spectrum stem 
from the BAE-CAP. Inside the island, they reach continuum accumulation points at 
the magnetic island O-point (MiO-CAP), while outside the island and far from the 
separatrix, they asymptotically approach the typical spatial dependence of the SAW 
continuum in a sheared magnetic field in the absence of the island. 

Understanding the modification of the SAW continuous spectrum due to the 
presence of a magnetic island has potential implications in explaining stability of Alfven 
instabilities in tokamak plasmas. Modes in the BAE frequency range have been observed 
in the Frascati Tokamak Upgrade (FTU) [27l [28] in the presence of an (m, n) = (-2,- 
1) magnetic island, where m and n are, respectively, the poloidal and toroidal mode 
number. A theoretical analysis showed that these modes can be interpreted as BAE 
modes, when thermal ion transit resonances and finite ion Larmor radius effects are 
accounted for, with good agreement of measured and calculated frequencies in the small 
magnetic island amplitude limit [29]. In fact, measured frequencies were found to depend 
on the magnetic island amplitude as well [27], consistently with the dependence of the 
BAE-CAP on the magnetic island size resulting from our theory. The modes were 
observed only when the magnetic island size was over a certain critical threshold [27] . 
Later on, similar observations have been reported in other tokamaks (see, for instance, 

X Ref. [25) analyzes this problem with a simple representation of the magnetic island as a straight 
flux-tube with non-circular cross-section. 
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the observations in HL-2A [3U]). This work derives and discusses the structures of 
the SAW continuous spectrum in tokamak geometry and in the presence of a finite 
size magnetic island. As mentioned above, the nonlinear modification of the BAE 
dispersion relation due to the presence of a magnetic island is not treated here and 
will be presented in a subsequent paper. The scheme of the paper is the following. In 
Sec. |2]the equilibrium magnetic field and the model equations are given; meanwhile, a 
convenient a coordinate system is defined inside and outside the magnetic island. In 
Sec. [3] we solve the problem for the continuous spectrum both numerically, in the whole 
region inside the magnetic island and analytically, near the O-point, obtaining the value 
of the MiO-CAP. Similarly, in Sec. H] we solve the problem for the region outside the 
magnetic island and compare the numerical solution with the asymptotic analytical 
solution far from the separatrix. The frequency of the BAE-CAP is also calculated as a 
function of the island size and mode numbers and parity. Sec. |5]is devoted to a summary 
of the novel obtained results and their possible application to the study of BAE in the 
presence of a magnetic island. 



2. Equilibrium and model equations 

2.1. Coordinate system 

We consider a tokamak geometry for a torus with major radius Rq. The equilibrium 
is made of an axisymmetric tokamak magnetic field with a component B tor in the 
toroidal direction and a component B po i in the poloidal direction 9t, plus an helical 
perturbation in the radial direction tt, generating a magnetic island. The magnetic 
island is located at a flux surface with minor radius r . Here, the subscripts T denote 
tokamak coordinates. We consider the region in the proximity of the rational surface 
of the magnetic island qr = % = rm s i/ n ish where m is i and n is i are, respectively, 
the poloidal and toroidal mode numbers of the magnetic island perturbation and 
q T = rTB tor / (RqB po i) is the safety factor. In this region, the toroidal magnetic field 
is assumed to be dominant and constant in space, B tor = B /j, and the poloidal 
magnetic field is considered to vary linearly with gj>. Here Bq is the axisymmetric 
magnetic field amplitude, 7 = i/l + £o/qi, and e = r /R . We adopt a slab model 
with coordinates (qT,u,(), where u is defined in Eq. [TJ by applying a rotation of the 
coordinates in the (6*t, Ct) plane, as shown in Fig. [TJ In this model, the coordinate 
C — (Ct + £q #t/<7o)/(<7o7 2 ) is the coordinate of translational symmetry. The equilibrium 
magnetic field is B = B^^tt + B u ^ p hU + B^^hC, where the axisymmetric field 
components are B uph = B e (qT — <?o)/(<?o7 2 ) an d -^Cp^ = -^0 (the physical components 



are defined in detail in Appendix A ). The constant-ip approximation is also adopted, 
assuming that the magnetic island field is B qT>p f l = B is i sinw, with B is[ constant. 

The flux surfaces of this equilibrium are labelled by ip, where the coordinates ip and 
u are defined by: 

V>= (gr-<?o)72 + M(cosM + l) , u = n isl {( T - q 6 T ) (1) 
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The X-points of the magnetic island are at (gy — qo,u) = (0,0) and (0, 2tt) and the 
O-point at (0,7r). M is a constant with value M = (qo\s\/n is i)(B is i/ B po ifl), determined 
by the condition • B = 0, where s is the magnetic shear and B po i is the poloidal 
magnetic field evaluated at the rational surface. 

In the slab model approximation, the plasma inside the magnetic island is a straight 
flux tube with length Z = jqoRo- The magnetic axis of the flux tube, directed along 
£, and the O-point of the magnetic island are at ip = and the separatrix is labeled by 
ip = i>sx — 2M. An appropriate set of cylinder-like coordinates (p, 9, Q is defined here 



to describe the region inside the magnetic island (see Appendix A ), with: 
P 



— — \[2a\) , = arccos {y/M{cosu + l)/i/>) (2) 

With these definitions, the magnetic axis is at p = and the separatrix radius is p = p sx , 
which corresponds to the magnetic island half-width, Wi S i, given by the Rutherford 
formula [3T] : 



W is i / B is i P. 



sx 



y qosn is iB po i fi r q^n is i 

with e defined below in this paragraph. The angle 9 is defined in the domain (0, 7r/2) and 
extended to (0, 7r) by reflection symmetry w.r.t. 9 = n/2. Further extension to (0, 2tt) 
is obtained by reflection symmetry for p -h- —p. With this definition, 9 has values 0, tt 
at the rational surface q = qo- We also point out that the flux surface's cross section 
in the (p, 9) plane and in the proximity of the O-point is an ellipse, with eccentricity 
e = 1 — Mn^ sl , y 2 j s 2 . Typical magnetic islands in tokamak experiments have values of 
eccentricity close to e ~ 1. Outside the island, the flux surfaces have the same topology 
of the original axisymmetric equilibrium, and therefore the coordinate system (p, u, £) 
is appropriate to describe the problem (see Eq. [Q. 
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2.2. Equilibrium magnetic field and model equations 

The equilibrium magnetic field in the coordinates (p, u, () is described by the contra- 



variant physical components: B p h = 0, B™ h = 2e B y/MP '/ '(<?o7 2 )> ^ph = Bo- Here P = 
a/L 2 + ( 1 — e) (sin 2 u) /4, where L = a/x 2 — (cos u + l)/2 and x = p/ p sx . The contra- 
variant physical components of a vector V, in a basis {<7 P , g u , g^} are defined here as the 
contra-variant components rescaled with the length of the correspondent basis vector, 
e -S- Vph = V P \dp\ (see |Appendix A| for further details). Similarly, in the coordinates 



defined inside the island we have B p h = 0, B d ph = 2e B \/My/ax/(q' 2 j 2 ), B ph = B . 
The function a is defined as a = sin 2 9 + cos 2 9(1 — e)F 2 , with F = y/1 — x 2 cos 2 9. 

The safety factor outside the island, q out , is defined as the flux-surface average of the 
ratio of the toroidal magnetic field and the sum of the poloidal and radial magnetic fields. 
The flux-surface average is performed separately on both sides of the magnetic island, 
with r > r and with r < r , along the coordinate y out = f° T r d9 T ^l + B 2 T>ph / 'Bj T ph . 
By definition, we have: 

r dy mt B tor 1 /- 2 ~ du 
q ° ut J2^ 5 p r 27m«y o (l-2v / MP/g ) 

Far from the separatrix (x > 1), the safety factor reaches asymptotically the linear 
behavior typical of the slab model without island. When the magnetic island amplitude 
M vanishes, this behavior is recovered for all x. In fact, in this case we have that 
2a/MP — > (qr — Qo) and therefore q out — > qr- An important result is that the safety 
factor at the separatrix [x — 1) is not qo, but there is a small difference, proportional 
to the magnetic island size. This difference is positive on one side of the island, where 
qr > qo, and negative on the side where qx < qo- The absolute value of this difference 
can be calculated at the separatrix from Eq. [31 in the limit M<1. We obtain: 



q ut, sx -Qo = - 1.27 VM (4) 

7T 

This result has important implications for the modification of the value of the BAE 
continuum accumulation point, which is located at the separatrix. This island- induced 
frequency shift of the BAE-CAP can be calculated in an approximated form as: 
Au ~ k\\VA = v A m(q outySX - q )/(q Z ) ~ 1.27y/M(v A /Zo){m/q ), where m is the 
poloidal mode number. In fact, we will see in Section 14.11 that the BAE continuum 
accumulation point has not the same frequency as without island, but a nonlinearly 
modified frequency which is proportional to \/M, consistently with the nonlinearly 
modified value of the safety factor given by Eq. HI 

Similarly, the safety factor inside the magnetic island can be defined as the average 
of Q = pB /(Z B e ) over 9, with B e = B e h Fy/l - e/(y/ap) (see also reference [25]): 



ph 

2 7 , x 2 1 , , 
7r|s|Vl-e ^VMn isl 

where we have introduced the complete elliptic integral of the first kind K(x) = 
J^ 2 d6/F. A similar definition for the safety factor was given in Ref. [26J, but our result 
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is a factor g smaller. The basic difference of our derivation from that of Ref. is 
that we take into account that the length of a magnetic island flux tube is Z ~ 27cq R , 
whereas in Ref. [26J a flux tube with length Zq — 2tcRq is considered, which is the major 
circumference of the tokamak. In other words, we point out that the magnetic island flux 
tube performs qo loops before closing on itself, and therefore all periodicity boundary 
conditions have to be imposed in its whole length, and not in the length of the tokamak 
circumference. This difference enters in particular the definition of the coordinate along 
the axis of the magnetic island flux tube (, and approximated in Ref. [26J by the toroidal 
coordinate Ct- 

In this work, we want to study the shear Alfven wave propagation near the resonant 
flux surfaces where the energy is absorbed by continuum damping; therefore, we focus 
on the dynamics of modes that are characterized by radial singular structures. The 
linear equation for radially localized shear Alfven modes in a compressible nonuniform 
tokamak plasma can be written in the form [22j [23] : 

^Vi^ + Z 2 V|, ViV|,0 - u2bae - cap VI<P = (6) 

where u A = va/Z . The frequency of the low-frequency SAW continuum accumulation 
point, delimiting the frequency gap of the Beta induced Alfven Eigenmode (BAE), is 
denoted as ujbae-cap- We adopt, here, the value of the BAE-CAP given in reference [21] 
and valid in the fluid limit, i.e. for frequencies much larger than the ion transit frequency: 



i mn t 

U BAE-CAP = -^\ — 7 + 7F ( 7 ) 



Ro V mi V4 Ti 

where Tj and T e are the ion and electron temperatures and is the ion mass. We 
focus on frequencies higher than ujbae-cap and consistently neglect kinetic effects 
associated with wave-particle resonances [24|. The operators Vn and Vj_ are parallel 
and perpendicular gradients with respect to the equilibrium magnetic field. 

In the following sections, the model equation for radially localized shear Alfven 
modes, Eq. [61 is written in the coordinate systems inside and outside the magnetic 
island and the numerical solution is provided and compared with analytical expressions. 

3. Solution inside the magnetic island 

3.1. Eigenvalue problem 

Here, we write the equation for SAW continuum modes, Eq. El in the coordinates inside 
the island, (x,0,£), introduced above. We consider continuum modes with the same 
helicity of the magnetic island: d/d( = 0. In this framework, the problem reduces to 
2 dimensions. The parallel and perpendicular differential operators have the following 
form: 

Z„V|| =VMn lsl F-^- 

Q2 

/4Vi = a— 
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0.2 0.4 0.6 0.8 1 



x 

Figure 2. Continuous spectrum Sl 2 (x) for the small eccentricity case, e < 1, 
corresponding to M ~ 1, plotted versus the radial position inside the island. Typical 
values of the equilibrium parameters have been chosen and n^i = 1. The O-point is at 
x = and the separatrix at x — 1. The branches of both even and odd eigenfunctions 
are shown. At the O-point, the continuum accumulation points (MiO-CAP) recover 
the value of the cylinder limit. At the separatrix, the first odd frequency branch 
tends to the original BAE-CAP (Q 2 = 0) and the others to the nonlinearly modified 
BAE-CAP (n 2 = Mn 2 sl ). 



and the boundary conditions of the problem are periodicity conditions in 9 on a 2n 
circle. By applying the differential operators in this explicit form, Eq. [6] can be written 
in the form of an eigenvalue problem: 

Q 2 F 8 „ 8 



J a F- 

\-Mn 2 sl a 89 89 



f = (8) 



where Q 2 = (u 2 — ^bae-cap)/ u a * s ^ e eigenvalue, and / = 8 2 4>/8 2 x is the 
eigenfunction. 

It is useful to write Eq. |8] in the form of a Schrodinger equation. The motivation is 
twofold. Firstly, we can use standard shooting method techniques to solve numerically 
the problem in the whole range inside the island: < x < 1. Secondly, we can solve 
analytically the problem in the proximity of the island O-point, approximating the 
potential, and compare the numerical solution with the analytical one near the O-point. 
The eigenvalue problem takes the form: 

8 2 £1 2 

h + T7ZT A inh - V m h = (9) 

isl 



89 2 Mn 2 



with A in = 1/F 2 , V in = (g" n - g' 2 n j \2g in )) / '(2g in ), g in = aF, and h = fy/gi^. Here the 
prime denotes the ^-derivative, e.g. g' = 8g/89. 

The solution of Eq. [9J can be found numerically with a shooting method code in 
the range < x < 1, for any value of the eccentricity < e < 1. The result is shown 
in Fig. [2] for e « 1, and in Fig. H] for e = 0.99. We choose typical values for the 
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equilibrium parameters, q = 2, s — 1, £ — 0.1, n is i = 1. The periodic boundary 
condition in 9 is satisfied by an infinite set of solutions, labeled by j = 1,2, .... Three 
main features describe the continuous spectrum inside the island. 1) The continuous 
spectrum branches Q 2 have continuum accumulation points at the O-point, named here 
MiO-CAP. 2) Even and odd eigenfunction have different eigenvalues, due to the non- 
uniformity of the magnetic field intensity along the field line. 3) At the separatrix, the 
continuum frequencies converge to two different BAE-CAP. The first odd eigenfunction 
has frequencies converging to the original BAE-CAP, whose value is the same as in an 
equilibrium without islands, correspondent to Q 2 = 0. All other branches converge to a 
nonlinearly modified BAE-CAP, whose value is Q = Mn 2 sl . 

The solution of Eq. [9] can also be found analytically near the O-point and near 
the separatrix, approximating the potential and obtaining the value of the continuum 
accumulation points (MiO-CAP and BAE-CAP). The analytical derivation of the MiO- 
CAP value is shown hereafter, and an estimation of the nonlinearly modified BAE-CAP 
is given in Section 14.11 



3.2. Analytic solution near the O-point 

The value of the continuum accumulation points at the O-point of the island can be 
calculated analytically from Eq. ((9]), by approximating the functions A in and Vi n near 
the O-point. In this limit, Eq. (Q becomes: 

where the potential is defined by: 

^n,oW = ef / 1 l ~ e)CQ f^ (2-ecos 2 g) - l) (11) 



^1 — ecos 2 i 

Two limiting cases are considered here: e C 1, corresponding to M ~ 1, and e ~ 1, 
corresponding to M <C 1. The former describes a magnetic island where the flux 
surfaces have circular cross section near the O-point. The latter is representative of 
typical size magnetic islands in tokamak plasmas. In the former case, the function 
a can be approximated by a ~ 1 and the magnetic field intensity is independent of 
9: B e ph ~ e \s\B p/(q rQj 2 ). This means that, in our model, the e C 1 case has a 
cylindrical symmetry around the O-point. In this case, the potential Vi n ^ vanishes and 
Eq. [10] has the following eigenvalues: 

^mo-CAP = Mn 2 slJ 2 = f/ql(0) (case e « 1) (12) 

where j is a natural number (j = 1,2...), and <2m(0) is the safety factor given in Eq. |5l 
calculated at x — 0. 

On the other hand, in the case of small island amplitude (e ~ 1), the potential 
can be written as Vi n> o — Vq + 5V, where Vq = — 1 and 5V = V\ H(9 e ff). Here 
V\ = 1/(1 — e) = s 2 j \Mn 2 sl ^ 2 ) and H(9 e ff) is a function with unitary value inside the 
set (|0| < 9 eff ) U(|tt -9\ < 6 ef f) and zero elsewhere, where 9 eff = (1/2) (y/M/\s\) 1/2 . 



2D continuous spectrum of SAW in the presence of a magnetic island 



10 



Since the potential V^ n> o is not constant in 6, in this case the eigenvalue Q 2 has a different 
value for even and odd eigenf unctions. However, an approximated value of Q 2 can be 
calculated by making some considerations on the potential shape. In fact, having 5V 
a support which is very localized in 9 for small amplitude islands, where odd modes 
are small, we can neglect its contribution and consider the problem determined by 
Vin,o = Vq = — 1. This gives the estimate for odd mode eigenvalues: 

^Iho-cap = Mn 2 sl (j 2 - 1) (case e ~ 1) (13) 
4. Solution outside the magnetic island 

In this section, we face the problem of SAW continuum modes, described by Eq. El 

outside the separatrix of a magnetic island, in the coordinate system, (x,u,Q- The 

parallel and perpendicular differential operators have the following form: 

i— d 
Z V|, =2VMn tsl L- 



du 

p2 Q2 



x 2 dx 2 

Only continuum modes with m = q^n are considered, where m and n are respectively 
the poloidal and toroidal mode numbers, in the tokamak coordinates 9t and (t- Those 
are the modes with the same helicity of the magnetic island: d/dC, = 0, and therefore, 
the problem reduces to 2 dimensions. By applying the differential operators in this 
explicit form, Eq. [6] can be written as in Sec. [3] as an eigenvalue problem: 
r tt 2 AL d „ d i „ 

1- + —Z-Z-LP 2 — f = 14 

VMn 2 isl P 2 du du\ J K J 

and, as for Eq. [8j this can be put in the form of a Schrodinger equation. The Schrodinger 
equation can be solved numerically for x > 1 with standard techniques and the numerical 
solution can be compared with the asymptotic limit, obtained analytically for x 3> 1. 
Moreover, analytical estimations of the frequency can be found near the separatrix. The 
eigenvalue problem takes the form: 
d 2 VL 2 

^n> h + iT7 2 A ™*h - V out h = (15) 
ou 2 Mnf sl 



with A out = 1/(4L 2 ), V out = {g" mi - g' 2 ut /{2g out ))/{2g out ), g out = 2LP 2 , and h = fy/g^. 

The periodic boundary condition in u is satisfied by an infinite set of solutions, 
labeled with mode number j (for go integer, j = n, with n the toroidal number in the 
tokamak coordinate). The solution of Eq. [15] is shown in Fig. [3] and Fig. HI for typical 
island size, M = 10~ 2 , and typical values of the equilibrium parameters, go = 2, s = 1, 
e = 0.1, n is i = 1. The asymptotic behavior at x ^> 1 is depicted in Fig. [3] for the case 
j = 1. Three main features describe the continuous spectrum outside the island. 1) 
The continuous spectrum for x 3> 1 reaches the asymptotical limit Q 2 ^, with a constant 
difference with the limit of the continuous spectrum in the limit of vanishing island, f^: 

Ql = n 2 < (q T - g ) 2 > u = Mn 2 sl3 2 (Ax 2 - 2) 

n 2 = n 2 lim (q T - g ) 2 = Mn 2 slJ 2 (Ax 2 ) 

M— >0 
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Figure 3. Continuous spectrum f2 2 , plotted versus x 2 , for M = 10~ 2 , corresponding 
to e ~ 0.99. Typical equilibrium parameters have been chosen, and n^i = 1. Outside 
the magnetic island, all the frequency branches converge to a nonlinearly modified 
BAE-CAP at the separatrix (x = 1), with value ^niBAE-CAP = ^ n lsi- On the 
other hand, for x 1, the asymptotic limit fi 2 ^ is approached. In our model, due 
to the constant--0 approximation, the solution without magnetic island fi 2 , is shifted 
with respect to the asymptotic limit: S7§ = fi 2 ^ + 2M. Even and odd modes have a 
negligible difference in frequency: Af2 2 ~ 0{M 2 ). 



This difference is due to the constant-?/' approximation which is made modeling the 
island magnetic field. A more realistic asymptotic behavior can be obtained by 
abandoning the constant-?/' approximation and considering an improved model where 
the island field decays with growing \qr — Qo\- 2) Even and odd eigenfunction have 
different eigenvalues, due to the non-uniformity of the magnetic field intensity along the 
field line, but the difference is negligible with respect to the absolute value. In the case 
depicted in Fig. [3J the difference between odd and even solutions is Afi 2 ~ 0(M 2 ). 
3) At the separatrix, all continuum frequencies converge to the nonlinearly modified 
BAE-CAP, whose value is proportional to the magnetic island half-width: 



This is consistent with the qualitative asymptotic behavior of the continuous spectrum 
near the separatrix, which can be calculated by knowing the value of the safety factor 
outside the magnetic island (see Eq. H]). Note that the nlBAE-CAP has the same 
frequency as the MiAE gap central frequency [25] . 

The nonlinear modification of the continuous spectrum and in particular the upward 
shift in frequency of the BAE-CAP has important implications in the study of the 
dynamics of discrete AE in tokamak equilibria in the presence of a magnetic island, 
as discussed in Sec. [5j The dynamics of discrete AE is not treated here, and will be 
discussed in a dedicated paper. The analytical treatment of the continuous spectrum 
problem near the separatrix is presented in next section. 




(16) 
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Figure 4. Continuous spectrum £l 2 (x), for M — 10~ 2 , corresponding to e ~ 0.99. 
Typical equilibrium parameters have been chosen, and 7ij S / = 1 . The region inside the 
island is at < x < 1, and the region outside the island at x > 1. The MiO-CAP 
are shown at the O-point (x — 0), and the linear and nonlinear BAE-CAP at the 
separatrix (x = 1). 



4-1. Analytical treatment near the separatrix 

Here, we focus on the region outside the magnetic island, in the proximity of the 
separatrix. We also consider the limit of e ~ 1, which is the case of a typical size 
magnetic island in tokamak plasmas. Under these assumptions, the function L and 
P can be approximated as L = P = y/1 — cosu/v2. Consequently, the eigenvalue 
equation, Eq. [151 takes the form: 

^h-W n2 h = (17) 

where the potential Wq2 is defined as: 

9 cos 2 u - (12 + m 2 ) cos u + (3 + 80 2 ) 



16(1 — cos it) 2 

and tt 2 = n 2 /(Mn 2 sl ). This equation can be studied as a Schrodinger equation with 
zero energy. In this framework, we look for values of Q 2 such that the solution satisfies 
the boundary conditions of periodicity in u on (0, 2tt). The potential Wn^ has positive 
second derivative for Q 2 < 3/4, and negative second derivative for Q 2 > 3/4. For 
Q 2 = 3/4, the potential has a constant value: W^p. = — (3/4) 2 . For the given boundary 
conditions, this means that solutions of the problem exist for Q 2 > 3/4 only. On the 
other hand, when Q 2 = 4/3 the potential has maximum value Wq2 = — 1 at u = it, and 
therefore does not admit a solution with mode j — 1. With these considerations, we can 
estimate that 3/4 < ^xbae-cap < 4/3 for the first mode number, j = 1. Numerically, 
we find that Q 2 ~ 1 for j = 1 and that the difference in frequency for other mode 
numbers and different parities is negligible. 
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5. Conclusions and discussion 

Understanding the structure of the SAW continuous spectrum is one crucial step in 
studying the stability properties of a tokamak plasma. In fact, one of the main damping 
mechanisms of SAW instabilities in non-uniform plasmas is continuum damping, which 
occurs near the magnetic flux surfaces where the frequency of the instability matches 
the SAW continuum frequency. In toroidal devices, a frequency gap in the continuous 
spectrum represents a range of frequency where discrete Alfven eigenmodes (AE) can 
grow unstable, practically unaffected by continuum damping. 

We have studied the radial structure of the SAW continuous spectrum in the 
presence of a magnetic island. Since magnetic islands are very common in fusion plasma 
experiments, for example in connection with tearing instabilities and sawtooth activity, 
the solution of this problem has important implications for understanding the stability 
regimes of AE in a fusion reactor. We have adopted a linear MHD model for radially 
localized SAW, where tokamak curvature is neglected except for the creation of the beta 
induced Alfven eigenmode (BAE) frequency gap, in the low-frequency range of the SAW 
continuous spectrum. We also assumed modes with the same helicity of the magnetic 
island. In this framework, the model reduces to 2 dimensions. The case of local singular 
structures with helicity different from that of the magnetic island is considered in a 
different paper [25]. Due to the time scale separation between island and the SAW 
dynamics, the island has been treated as a static helical distortion of the equilibrium. 

We have shown that the SAW continuous spectrum of a finite-/? tokamak plasma 
is modified by the presence of a magnetic island. In particular, the BAE continuum 
accumulation point (BAE-CAP) is shifted in space from the island rational surface, to 
the separatrix flux surface position. Moreover, we have shown that the nonuniformity of 
the magnetic field intensity along the field lines generates a frequency splitting between 
modes with even and odd parity. This splitting frequency has a finite value inside the 
island and is negligible outside for typical island sizes. Inside the magnetic island, we 
have found that the continuum frequency branches converge at the O-point to several 
magnetic island induced CAP (MiO-CAP). At the separatrix, the frequency of the 
BAE-CAP is split in two values. The continuum branch corresponding to the first odd 
eigenfunction inside the island converges to the original BAE-CAP with same frequency 
as without magnetic island; and the other branches converge to a nonlinearly modified 
BAE-CAP, with a higher frequency given by Eq. [16j Outside the magnetic island, 
all continuum branches at the separatrix converge to the nonlinearly modified BAE- 
CAP. Far from the separatrix and outside the magnetic island, the continuum frequency 
converges asymptotically to the value calculated in an equilibrium with no island. 

Our results have two main implications in the study of the dynamics of Alfven 
eigenmodes (AE). Inside the island, new magnetic island induced AE (MiAE) could 
exist as bound states, essentially free of continuum damping, provided that plasma 
equilibrium effects and free energy sources can drive and bind them locally (see also 
reference [25]). Outside the island, BAE are expected to have an eigenfunction radially 
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peaked in the proximity of the magnetic island separatrix, namely at the BAE-CAP, 
other than at the rational surface. The BAE frequency is expected to have higher 
values when a magnetic island is present, consistently with the BAE-CAP frequency 
modification described by Eq. [16j Moreover, while inside the magnetic island the 
pressure and temperature profiles are flattened, outside the magnetic island, in the 
proximity of the separatrix, the equilibrium pressure gradients are highly increased. 
This may provide a free-energy source for BAE, which are expected to grow unstable in 
this configuration, even in absence of fast particle drive. BAE are generally damped by 
thermal ion Landau damping; therefore, we may expect to detect BAE when the value 
of the pressure gradient is over a certain threshold, namely when the magnetic island 
exceeds a certain size. The dispersion relation and the drive and damping mechanisms 
of BAE in the presence of a magnetic island will be treated in a different paper. 

Experimental observations of Alfvenic modes in the presence of a magnetic island 
have been recently reported in several tokamaks devices (see for instance the observations 
in FTU [271 128] and in HL-2A [30J ) , even in purely Ohmic heated plasmas. The 
fluctuations frequency range confirms that they are BAE nonlinearly interacting with 
the magnetic island, characterized by a frequency shift dependence on the island 
size consistent with that of the nonlinearly modified BAE-CAP. Moreover, BAE were 
observed only above a certain threshold in the magnetic island amplitude. 

The dependence on the magnetic island size of the continuous spectrum and in 
particular the MiO-CAP and the BAE-CAP suggests the possibility of using the MiAE 
and BAE frequency scalings as novel magnetic island diagnostics in a fashion similar 
to other commonly used Alfven spectroscopy techniques J32J [331 [31] . The radial MiAE 
localization at the center of the island makes them more difficult to detect by external 
measurements than BAE and makes the use of internal fluctuation diagnostics, such as 
electron cyclotron emission (ECE) and soft X-Rays, necessary. On the other hand, BAE 
are easier to detect with external diagnostic techniques, and therefore, the possibility of 
using BAE frequency scalings as a magnetic island diagnostic can be considered more 
feasible. 
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Appendix A. Coordinate metric 

Here, we provide the metric for the coordinates used in this work outside the magnetic 
island (p, u, () and inside the magnetic island (p, 8, (). This is necessary to calculate the 
differential operators in the equation for the dynamics of continuum modes. 

The domains of the coordinates (p, u, C), describing the region outside the island, 
are p sx = W is i < p < oo, < u < 2rc and < ( < 2n, where ( is the coordinate of 
translational symmetry for both equilibrium and perturbations, and periodicity in u is 
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assumed for the perturbations. The gradients of these coordinates are: 
Vp = L/xq T — y/1 — esinu/(2x) u = (P/x) p 
Vm = u/po 

vc = C/z 

where p — r o/{ < lo n isil)i Zq = 7^0-^0 an d we use the notation V = V /V for a versor 
(unit vector). Here P = \/ L 2 + (1 — e)(sin 2 -u)/4, where L = x 2 — (cos-u + l)/2 and 

X = p/psx- 

A covariant basis (g p , g u , g^) can also be defined as orthogonal to the contra- variant 
basis (g p = Vp, g u = V«,g^ = VC), namely satisfying: g l gj = <5}. With this notation, 
any vector can be written as V = ^ V l <7z = £V Vph9ii where we call V* h the physical 
components of the vector, V£ h = V l \gi\. The Jacobian for the coordinates outside the 
island, defined as the determinant of the metric tensor g^ = g t ■ gj, is g out = p^Z^x 2 / L 2 . 
With these definitions, the laplacian of a function / for "radially" localized modes, can 
be written in the coordinates (p, u, () as: 

y/g^Tt dp\ v dp J pf x x 2 dx 1 

where we have used g pp = g p ■ g p = P 2 /x 2 . 

The domains of the coordinates (p, 9), for the region inside the island, are < p < 
Psx = W is i, < 9 < 2rc and < ( < 2tc, where ( is the coordinate of translational 
symmetry as above, and periodicity in 9 and ( is assumed for the perturbations. The 
gradients of the coordinates p and 9 are: 



Vp = sin 9 q T — cos 9y/l — eF u = \fa p 

V# = (cos 9q T + sin9Vl-eFu)/p = VbO/p 

Moreover we have the relation Vp ■ V6* = c/p. Here F = y/l — x 2 cos 2 9, a = 
sin 2 9 + cos 2 #(l — e)F 2 , and c = cos#sin#(l — (1 — e)F 2 ). The Jacobian for the 
coordinates inside the island, is g in = p 2 x ZqX 2 /(F 2 (l — e)). Using these definitions, and 
the fact that g pp = g p ■ g p = a, we write the laplacian for "radially" localized modes in 
the coordinates (p, 9, () as: 

V 2 f ~ —^L (A.2) 



